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We investigate exact solution in higher dimensional Husain model for a null fluid source with 
' pressure p and density p are related by the following relations (i) p — kp, (ii) p = kp— (variable 

, modified Chaplygin) and (iii) p = kp" (polytropic) . We have studied the nature of singularity in 

. gravitational collapse for the above equations of state and also for different choices of the of the 

(N . 

parameters k and B namely, (i) k = 0, B — constant (generalized Chaplygin), (ii) B — constant 

^ l' (modified Chaplygin). It is found that the nature of singularity is independent of these choices of 

U ■ different equation of state except for variable Chaplygin model. Choices of various parameters are 

' shown in tabular form. Finally, matching of Szekeres model with exterior Ifusain space-time is done. 
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CJ ; I. INTRODUCTION 



Gravitational collapse is one of the most important problem in classical general relativity. Usually, the 
formation of compact stellar objects such as white dwarf and neutron star are preceded by a period of collapse. 
The study of gravitational collapse was started by Oppenheimer and Snyder [1] . They studied collapse of dust 
' with a static Schwarzschild exterior while interior space-time is represented by Friedman like solution. Since 
, then several authors have extended the above study of collapse of which important and realistic generalizations 
' are the following: (i) the static exterior was studied by Misner and Sharp [2] for a perfect fluid in the interior, 
ly-^ I (ii) using the idea of outgoing radiation of the collapsing body by Vaidya [3], Santos and collaborations 
■ [4-9] included the dissipation in the source by allowing radial heat flow (while the body undergoes radiating 
collapse). Ghosh and Deskar [10] have considered collapse of a radiating star with a plane symmetric boundary 
(which has a close resemblance with spherical symmetry [11]) and have concluded with some general remarks. 
' Also Cissoko et al [12] have studied junction conditions between static and non-static space-times for analyzing 
, gravitational collapse in the presence of cosmological constant. 



So far most of the studies have considered in a star whose interior geometry is spherical. But in the real 
astrophysical situation the geometry of the interior of a star may not be exactly spherical, rather quasi-spherical 
5h \ in form. Recently, solutions for arbitrary dimensional Szekeres' model with perfect fluid (or dust) [13] has 
been found for quasi-spherical or quasi-cylindrical symmetry of the space-time. Also a detailed analysis of 
the gravitational collapse [14, 15] has been done for quasi-spherical symmetry of the Szekeres' model. It has 
also been studied junction conditions between quasi-spherical interior geometry of radiating star and exterior 
Vaidya metric [16]. 

In 1996, Husain [17] gave non-static spherically symmetric solutions of the Einstein equations for a null 
fluid source with pressure p and density p related by p = kp. Wang et al [18] has generalized the Vaidya 
solution which include most of the known solutions to the Einstein equation such as anti-de-Sitter charged 
Vaidya solution. Husain solution has been used to study the formation of a black hole with short hair [19] and 
can be considered as a generalization of Vaidya solution [18]. Recently, Patil et al [20, 21] have studied the 
gravitational collapse of the Husain solution in four and five dimensional space-times. 

Recent observational data shows anisotropy in the Cosmic Microwave Background Radiation (CMBR) [22] 
and the evidences from type la Supernovae SN 1997H redshift [23 - 25] suggest that the universe is flat and 
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is undergoing at present an accelerating phase preceded by a period of deceleration. But known ordinary 

(baryonic) matter or radiation can not explain these observational facts, there should be a significant amount 
of extraordinary non-baryonic matter (dark matter) and energy (dark energy). In last few years, several 
models have been suggested to incorporate the recent observational evidences and among them, a single 
component perfect fluid having exotic equation of state, known as Chaplygin gas [26, 27] is of great interest. 
A further generalized form is known as modified Chaplygin gas (sec cq. (16) with B — constant with equation 
of state) (MCG) [28]. The MCG model can interpolate states between standard fluids at high density and 
constant -ve pressure at low density. For 7 = 1/3, the MCG model describes evolution from radiation epoch 
to ACDM era at late time (where the fluid has constant energy density and behaves as a cosmological constant). 

Further, from phenomenological point of view MCG model is interesting and can be motivated by the brane 
world interpretation [29]. Moreover, this model is consistent with various classes of cosmological tests namely 
gravitational lensing [30], gamma-ray bursts [31] as well as the above mentioned observations. Also for low 
energy density, the equation of state is that of a polytropic gas [32] with -ve index. Therefore it is possible 
to have astrophysical implications of the present model with an alternative way of restricting the parameters [33]. 

Moreover, the Chaplygin gas model produces oscillations or exponential blow up of the matter power 
spectrum that are inconsistent with observation [33]. Variable Chaplygin gas model interpolates between 
a universe dominated by dust and a quiessence dominated universe described by the constant equation of 
state. Also the model corresponds to a Born-Infeld tachyon action [29]. Recently, the model parameters were 
constrained using the location of peaks of the CMB spectrum [32] and SNIa data [34]. Also the variable 
Chaplygin gas model is constrained using SNela gold data [35]. 

The study of higher dimensional space-time model can be justified from two perspective. Firstly, inspiration 
comes from supcrstring theory and supcrgravity theory where one must have higher dimensional space-time. 
Also the idea of brane world scenario needs extra dimension for bulk space. Secondly, in recent past it has 
been shown that gravitational collapse in higher dimension has some peculiarity at the fifth dimension. So for 
the present model we have investigated collapse dynamics in higher dimensional model. 

In this paper, we give the generalization of Husain solution in (n -|- 2)-dimensional spherically symmetric 
space-time for a null fluid source with pressure p and density p are related by the following relations (i) 
p = kp, (ii) p = kp— (variable modified Chaplygin) and (iii) p = kp" (polytropic) [36]. For A; = 0, (ii) 

reduces to variable Chaplygin model [35]. We have also considered the particular cases namely, (i) k = 0, B = 
constant (generalized Chaplygin) [27], (ii) B = constant (modified Chaplygin) [28]. We have considered the 
interior space-time by Szekeres' model while for exterior geometry we have considered generalization of Husain 
space-time and is presented in the appendix. 



II. BASIC EQUATIONS FOR HUSAIN SOLUTION IN (n + 2)-DIMENSIONAL 
SPHERICALLY SYMMETRIC SPACE-TIME 

The metric ansatz in (n -|- 2)-dimensional spherically symmetric space-time can be taken as 

ds^ = - (1 - ^^J:^) dv^ + '^dvdr + r^d^l (1) 

where as usual r is the radial co-ordinate (0 < r < 00), the null co-ordinate v {—00 < v < 00) stands for 
advanced Eddington time co-ordinate such that r decreases towards the future direction along any v =constant 
trajectory, m{v, r) gives the gravitational mass inside the sphere of radius r and dfl'^ is the line element on a 
unit n-sphere. 

For the matter field, we have two non-interacting components namely the Vaidya null radiation and a perfect 
fluid having form of the energy momentum tensor 

rp _ rp{n) I rp{m) ( C,\ 
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with 

T^u^ = <^l^lu (3) 

and 

T^u^ = iP + pW^Vu + luV^) + P9^. (4) 

In the comoving co-ordinates {v,r, 61,62, the two eigen vectors of energy- momentum tensor namely 1^ 
and riv are linearly independent future pointing light-like vectors (null vectors) having components 

or explicitly, 



= (1, 0, 0, 0) and 77^ = Q (l - ,-1,0, 0^ 



(5) 



and they satisfy the relations 



Ixl^ = VxV^ = 0, hv^ = -1 (6) 



Here p and p are the energy density and thermodynamic pressure while a is the energy density corresponding 
to Vaidya null radiation. For p = p = 0, the matter distribution corresponds to null radiation flowing in the 
radial direction. When a = 0, then matter reduces to degenerate type I fluid [37] (In type I fluid, the energy- 
momentum tensor has a time-like eigen vectors and 3-space-like eigen vectors). In order to satisfy the energy 
conditions by the stress-energy tensor (2), the restrictions on p, p and a are 

(a) Weak and strong energy conditions: cr > 0, p > 0, p > (7) 

(b) Dominant energy conditions: a > 0, p > 0, p > p (8) 
The non-vanishing components of the Einstein field equations (choosing 8ttG = 1 = c) 

for the metric (1) with matter field having stress-energy tensor given by (2) are 

nm' m" , nm 

P=^^ P = -^i ^'^da=— (9) 

where an overdot and dash stand for differentiation with respect to v and r respectively. In order to satisfy 
the energy conditions (7) and (8) we have the following restrictions on m from the field equations (9): 

(i) m' > 0, to" < and (ii) m > 0. 

The condition (i) means that the mass function either increases with r or is a constant - a natural physical 
requirement. The second restriction (ii) implies that matter within radius r increases with time - an implosion. 

We shall now solve the above field equations (9) for the following equation of state: 

(a) Linear equation of state, 

(b) Variable modified Chaplygin gas model. 
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(c) Polytropic fluid. 

Case (a): Linear equation of state: 

Here the linear equation of state for the perfect fluid is taken as 

p = kp, {k, a constant) 



(10) 



Due to dominant energy condition k should be less than unity, but for other energy conditions k(> 0) is 
unrestricted and may even be greater than unity. Now using (10) in the Einstein field equations, we have the 
differential equation in m, 



m" = - 



knm' 



(11) 



for which the explicit solution for m is 



m{v,r) 



f{v)+g{v)logr, nk = 1 



(12) 



where the integration functions f{v) and g{v) are not totally arbitrary but are restricted by the energy 
conditions. 

The expression for physical quantities are (for nk ^ 1) 



p = kp 



nkg{v) 

2^n(fe+l) 



(13) 



(7 = 



f{v) 



9{v) 



{nk - l)r"fc-i 



(14) 



From equation (13) we must have g{v) > for positive pressure and density. The dominant energy condition 
leads, for nfc ^ 1, to 



m = f{v) 



{nk — l)r' 



,nk—l 



> 



Physically this means that the matter within a radius r increases with time, which corresponds to an implosion. 
This condition is satisfied if f{v) > and { either (a) g{v) > 0, fc < 1/n or (b) g{v) <0,k > 1/n}. 

Therefore the solution of the Einstein equations with matter field given by equation (2) can be written as 



fjv) 



+ 



9{v) 



{nk - l)r^{k+i)-2 



dv^ + 2dvdr + r'^dnl 



(15) 



This is termed as Husain solution in higher dimension. Similar to four dimensional Husain metric, g{v) = 
corresponds to usual higher dimensional Vaidya metric while for fc = — l), the solution represents Charged 
vaidya metric. Further, one should note that the above metric (15) is asymptotically flat for k > 1/n and is 
cosmological for A; < 1 /n. 



On the otherhand, for the solution with k = 1/n, the energy conditions are not always satisfied for all r 
as m = f{v)+g{v) log r becomes negative (i.e., a <0) for very small r. Hence we shall not discuss this solution. 
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Case (b): Variable modified Chaplygin gas: 

Here the mathematical relation (equation of state) between pressure and density is chosen as 

P = kp-^ (16) 

with fc(> 0), a{> 0) being constants and B{v), a function of v alone. Using (16) in the field equations (9), 
the differential equation in m becomes 

„ _ nkm' B{v){n - i)«+V"(«+i)-i 
^ ~ r ^ (nm')" 

which has the solution (for a ^ 1, nk ^ 1) 

, X {9{v)}^ 1-nk 1 1-nk 

m{v,r) = f{v)- ^F,[ --—,1 



(nA;-l)r"'=-i ^n(l + fc)(l + Q!)' 1 + a n(l + fc)(l + a)' 

{l + k)n^+°'g{v) ^ 



(18) 



Here 2-F1 is the usual hypergeometric function and f{v) and g{v) are as before arbitrary integration functions. 
If one chooses B to be constant then the above equation of state (16) corresponds to modified Chaplygin gas 

while k = 0, B =constant represents generalized Chaplygin gas and k = 0, B ^constant, a = 1 corresponds to 
pure Chaplygin gas. Further, for B = 0, one gets back the solution (15) for linear equation of state. 

Ccise (c): Polytropic fluid: 

The form of equation of state for such fluid is 

p = kp", {k, a constants) (19) 
So from the Einstein equations the differential equation becomes 

m" = - (20) 

^^a — lj,na—n+l ^ ' 

As before, the solution in terms of hypergeometric function can be written as (a ^ 1) 

m{v,r)=!{v)+r{g{y)')r^ ^F^\— — ,1 + — j^— ^ ] (21) 

n(l — a) 1 — a n(l — a) 9\v) 

This solution can be obtained from the solution (18) as a particular case fc = 0, B{y) — > — fc, a. — > —a. 



III. STUDY OF GEODESIC AND EXISTENCE OF NAKED SINGULARITY IN HUSAIN 

SPACE-TIME 

We shall discuss the existence of naked singularity in Husain space-time by studying radial null geodesies. In 
fact, we shall examine whether it is possible to have outgoing radial null geodesies which were terminated in 
the past at the central singularity r = 0. The nature of the singularity (naked singularity or black hole) can be 
characterized by the existence of radial null geodesies emerging from the singularity. The singularity is at least 
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locally naked if there exist such geodesies and if no such geodesies exist it is a black hole. 

The equation for outgoing radial null geodesies can be obtained from equation (15) (for linear equation of 
state) by putting ds^ = (with dfl"^ = 0) as 

2^ = 1 _ M + 9{v) .221 

dv r"-i {nk - l)r»(fc+i)-2 ^ 

It can be seen easily that r = 0, v = corresponds to a singularity of the above differential equation. 
Suppose X = ^ then we shall study the limiting behaviour of the function X as we approach the singularity at 
r = 0, V = along the radial null geodesic identified above. If we denote the limiting value by Xq then 



Xo =lim X = Urn ^ = lim 3? = 2 



where the arbitrary functions f{v) and g{v) are chosen as 



1 A yn-l , fi, yn{k+l)-2 ^ 



(23) 



f{v) = Aw"-i 
g{v) = ^v"(fc+i)-2 



(24) 



Note that only the above choices of / and g have given non-zero contribution to the limiting process given 
in equation (23). Otherwise the tangent to the radial null geodesic will have trivial solution. 

If the above polynomial in Xo has at least one positive real root then it is possible to have a radial null 
geodesic outgoing from the central singularity. More explicitely, Xq is a root of the following equation in X: 



{nk - 1) 



Xn(k+r,-i _ ^x-2 = (25) 



As exact analytic solution for X is not possible for general n and k so we consider the following particular 
case, namely n = 2, k = (i.e., dust solution in four dimension). The above root equation simplifies to a 
quadratic equation 



XX^ + (M - 1)X + 2 = i.e., Xo = ^-^^^^il^^-^r- ^ (26) 

Hence, if /x < 1 and 1 — /x > 2\/2A then Xq is always positive. As f{v) and g{v) are arbitrary function so 

we can choose A and /i appropriately so that Xq is positive definite. Hence, dust collapse in four dimension 
always leads to naked singularity. On the other hand for any n{> 2) we study the roots by numerical methods. 
Table I shows the possibility of radial null geodesic in various dimensions for different choices of parameters 

involved. We see that for fc = 0, there do not exist any radial null geodesic above 5D for various choices of 
other parameters, while for k = 1/3 or 2/3 it is possible to have radial null geodesies in 4D and higher dimensions. 

On the otherhand, if we use modified Chaplygin gas instead oip = kp then also the above root equation (25) 
remains unchanged provided 

f{v) = At;"-i 

^(.y) = ^l+a^(l+a){n(l+fe)-2} 

is chosen. Thus the nature of singularity of Husain model is independent of the choice of the matter in the form 
of perfect fluid having equation of state (i) p = kp, {u)p = — (pure Chaplygin), (iii) P = (generalized 

Chaplygin), (iv) p — kp ~ (modified Chaplygin) and (v) p ~ kp"' (polytropic). However, if we consider 
variable modified Chaplygin gas having equation of state p = kp — then the root equation (25) modifies to 

M ^„(i+fc)-i ^-nk 1 ^ , l-nk 



(nk-l) ^n(l + A;)(l + a)' 1 + a' n{l + k){l + a) 
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(1 + fc)(n/x)i+«X"(i+'=)(i+«) ^ 



X 



(27) 



where B(v) = Bov~'^^^^"\ Table II shows the nature of the singularities in presence of variable modified 
Chaplygin gas. Clearly, the nature of the roots will be different and hence variable modified Chaplygin gas has 
an effect on the collapsing process of the Husain model. For the choice a = 1, we have peculiar behaviour in 
5D and 6D, where there exist no radial null geodesic but it is possible to have radial null geodesic in 4D and 
higher dimensions. Also for k = 0, we have similar behaviour as in table I. 



A 




k 




Positive roots (Xq) 












4D 


5D 6D 


7D 


8D 


lOD 


0.1 


0.1 


2/3 


1.47479 


1.67881 1.67375 


1.62512 


1.56933 


1.47331 


0.01 


0.001 


1/3 


2.05205, 
91.7181 


2.12534, 

4.96352 
36.8271 


17.1871 


9.99197 


5.40875 


0.01 


0.1 


1/3 




1.246 


1.31917 


1.32766 


1.30382 


0.1 


0.001 





2.82809 
7.07191 










0.01 


0.1 





2.27998 
87.72 










0.001 


0.1 





2.22774 
897.772 


2.81434 
6.50088 









Table I: Nature of the roots (Xq) of the equation (25) for various values of parameters involved. 



A 




k 


a 




Positive roots (Xq) 














4D 


5D 


6D 7D 


8D 


lOD 


0.01 


0.001 


2/3 


1 


2.08069 


2.11867 


2.2056 3.17953 


2.95017 


2.36549 


0.1 


0.01 


1/4 


1 


3.08023, 
6.33532 




39.0619 


25 


10 


0.1 


0.01 


1/3 


1 


3.40688 
5.31946 




17.2107 


9.9992 


5.40875 


0.01 


0.01 





1 


2.925 
7.05688 










0.1 


0.01 





0.5 


2.93501 
7.04588 










0.01 


0.001 





0.5 


2.12406 
97.856 


2.18144 
8.72369 


2.40461 
3.38481 







Table II: Nature of the roots {Xq) of the equation (27) in variable modified Chaplygin gas for various values 
of parameters. 
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IV. CONCLUSION 



We have obtained exact solution in higher dimensional Husain model for null fluid source with equation of 
state between p and p are chosen as (i) p = kp, (ii) p = kp — (variable modified Chaplygin) and (iii) 
p = kp°' (poly tropic). We have seen that solutions for one equation of state can be obtained from the other by 
suitable transformation of the parameters. An interesting result has been obtained in studying the gravitational 
collapse of higher dimensional Husain space-time — the final state of collapse is independent of the choice of 
equation of states namely, linear, generalized Chaplygin, modified Chaplygin, polytopic. However, the result 
will change only if variable modified Chaplygin gas is considered {B is a function of v). From the table we see 
that in higher dimension, the value of a becomes insignificant to determine the positive root Xq. Finally, we 
have investigate the matching conditions of Husain model with interior Szekeres space-time in the appendix. 



APPENDIX 

Junction Conditions: Matching between Quasi-Spherical Szekeres model and Husain Space- 
Time 

Let E be a time- like (n + l)-dimensional hypersurface which divides (n -|- 2)-dimensional space-time into two 
distinct (n -|- 2)-dimensional manifolds V~ and The modified version of Israel [38] by Santos [4, 5] will 
be used to find the junction conditions. The geometry of the interior manifold V~ is given by the Szekeres 
space-time 

n 

dsl = -dt" + e^'^dr^ + e^^ ^ dx\ (28) 

i=l 

where a and /3 arc functions of all the (n + 2) space-time variables. The stress-energy tensor of a non-viscous 
heat conducting fluid has the expression [39] 

Tf^v = (p + V)'^p^v + V9ixv + (lixU^ + (IvU^ (29) 

where p, p, are the fluid density, isotropic pressure and heat flow vector. We take the heat flow vector to be 
orthogonal to the velocity vector i.e., q^u^ = 0. For comoving co-ordinate system we choose u'^ = (1, 0, 0, 0, 0) 
and qi^ = {0,q,qi,q2, ...,qn) where q = q{t,r,xi, ...,Xn) and qi = qi{t,r,x\, ...Xn),i = 1,2, ...,n. Now from the 
non-vanishing components of the Einstein field equation 

we have the solutions of (3 and a as [39] 



e'' = ii(i,r)e'^('''^^'^=--^") (30) 

and 

D{t,r) ^^^> 
where R and D are functions of t, r only. Here R and D satisfy 



2RR +(n- -D^) + - pR"" = (n- l)/(r) (32) 

n 

and 

Rb = f{r) e-2« (33) 

where /(r) is the arbitrary function of r. 



The function v satisfies the equation 
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e-'' = ^(r)^a;2+^B,(r)xi+C(r) (34) 



where A, Bi, C are arbitrary functions of r alone with the restriction 



- A AC = /(r) - 1 (35) 



Now from Einstein's field equations we have the components of heat flux as 



and 



77 

9=-De-« (36) 



qi = -^a,,e-^ (37) 



However, from equation (33) wc note that as R and D are functions of t and r only, so a is independent of 
the space co-ordinates x,; 's {i = 1,2, ...,n) i.e., a^^ = 0, V i = 1,2, ...,n. Hence from equation (37) we have 
Qi = and from equation (36) we have seen that q = q{t, r) i.e., g is a function of t and r only. Thus only radial 
heat flow is possible for the above choice of the metric as we consider. 

However, for the co-ordinate transformation [13] {xi,X2, •••,a;„) — > (^1,^2, ■■■,6n)'- 

xi = sin9nsin9n-i sin62Cot^6\ 

X2 = cos6nSin6n-i sin62cot^0i 

X3 = cos0n-isin9n-2 sin92Cot^6i 

(38) 



Xn-i = cos6zsin62cot\6\ 

Xn = COs92COt^9i 



the metric in V becomes 

dsl = -de + e""dr2 + i e^Pcosec\ei/2) d^l (39) 
where dVt^ is the metric on unit n-sphere. 

For the exterior space-time we consider (n -|- 2)-dimensional Husain space-time having metric 



dsl 



_ m{v,z) 



dv'^ + 2dvdz + z'^dQl (40) 



where m{v,z) is given in equation (12). 

The intrinsic metric on the boundary S is given by 



dsl = -dT^ + A^{T)dnl 



(41) 
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Now Israel's junction conditions (as described by Santos) are 

(i) The continuity of the Une element i.e., 

where ( )e means the value of ( ) on E. 

(ii) The continuity of extrinsic curvature over S gives 

[K,^] = K±-Kr.=o, 
where the explicit form of extrinsic curvature has the expression 



(42) 



(43) 



(44) 



Here = (r, a;i, a;2, a;„) are the intrinsic co-ordinates toE, xj, fT = 0,l,2,...,n-|-l are the co-ordinates 
in and are the components of the normal vector to S in the co-ordinates x± • 

Now for the interior space-time described by Szekeres metric the boundary of the interior matter distribution 
(i.e., the surface S) will be characterized by 

f{r,t)=r-r^=0 (45) 
where is a constant. As the vector with components is orthogonal to S so we take 

n; = (0,e«,0,...,0). 

So comparing the metric ansatzes given by equations (28) and (41) for dr = we have from the continuity 
relation (42) 



= 1, Air) = i e^cosec^ (6*1/2) on r = r„ 
dr 2 

Also the components of the extrinsic curvature for the interior space-time are 



(46) 



K^^ = and K^^g^ = cosec^eiKg^g^ 



J/3'e2^-"cosec^(^i/2) 



(47) 



On the other hand for the exterior Husain metric described by the equation (40) with its exterior boundary, 
^iven by 



f{z,v) =z- z^{v) = 
So the unit normal vector to S is given by 



(48) 



m(v,z) d,z\ 
1 -2— 1 



yU—l 



dv J 



dz 
dv 



,1,0,...,0 



(49) 



and the non-zero components of the extrinsic curvatures are 



l)mv 



22;" 



+ 



V dm 
2^"-i ~dz 



and 



/) = cose(?9iKt (, = .... = 



-zz + zv 



V [I q- 



(50) 



(51) 
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where 

. _2 -, m dz , . 

V =1 — 2— and = tA 

z^-^ dv ^ 

results from the contimiity of metric ansatz (40) on E (here an overdot stands for differentiation w.r.t. r). Now 
the junction condition due to the continuity of the extrinsic curvature results 



m{v,z) = ( - e^cosec^ (6ii/2) 



1 + icosec4(0i/2) (e2^/j2 - e^/^-^"/?'^) 



(52) 



and 



On the boundary vanishing of the isotropic pressure does not imply the vanishing of the heat flux. Thus for 
a quasi- spherical shearing distribution of a collapsing fluid, undergoing dissipation in the form of heat flow, 

the isotropic pressure on the surface of discontinuity E does not balance the radiation. Hence in the absence 
of isotropic pressure there may still be radiation on the boundary and the exterior space-time will still be 
Husain space-time. 



Acknowledgement : 



One of the authors (UD) is thankful to the authority of Institute of Mathematical Sciences, Chennai, India 
for providing Associateship Programme under which part of the work was carried out. Also UD is thankful to 
CSIR, Govt, of India for providing research project grant (No. 25(0153)/06/EMR-II). 



References: 



[I] J. R. Oppenhiemcr and H. Snyder, Phys. Rev. 56 455 (1939). 
[2] C. W. Misner and D. Sharp, Phys. Rev. 136 b571 (1964). 

[3] P. C. Vaidya, Proc. Indian Acad. Sci. A 33 264 (1951). 

[4] N. O. Santos, Phys. Lett. A 106 296 (1984). 

[5] N. O. Santos, Mon. Not. R. Astr. Soc. 216 403 (1985). 

[6] A. K. G. de Oliveira, N. O. Santos and C. A. Kolassis, Mon. Not. R. Astr. Soc. 216 1001 (1985). 
[7] A. K. G. de Oliveira, J. A. de F. Pacheco and N. O. Santos, Mon. Not. R. Astr. Soc. 220 405 (1986). 
[8] A. K. G. de Oliveira and N. O. Santos, Astrophys. J. 312 640 (1987). 

[9] A. K. G. de Oliveira, C. A. Kolassis and N. O. Santos, Mon. Not. R. Astr. Soc. 231 1011 (1988). 
[10] S. G. Ghosh and D. W. Deshkar, Int. J. Mod. Phys. D 12 317 (2003). 

[II] S. G. Ghosh and D. W. Deshkar, Gravitation and Cosmology 6 1 (2000). 

[12] M. Cissoko, J. Fabris, J. Gariel, G. L. Dcnmat and N. O. Santos, 3r-gc/9809057; S. M. C. V. Gonalves, 
Class. Quantum Grav. 18 4517-4530 (2001). 

[13] S. Chakraborty and U. Debnath, Int. J. Mod. Phys. D 13 1085 (2004); (/r-gc/0304072. 
[14] P. Szekeres, Phys. Rev. D 12 2941 (1975). 

[15] U. Debnath, S. Chakraborty and J. D. Barrow, Gen. Rel. Grav., 36 231 (2004); (?r-gc/0305075. 

[16] U. Debnath, S. Nath and S. Chakraborty, Gen. Rel. Grav. 37 215 (2005 ). 

[17] V. Husain, Phys. Rev. D 53 R1759 (1996). 

[18] A. Wang and Y. Wu, Gen. Rel. Grav. 31 107 (1999). 

[19] J. D. Brown and V. Husain, Int. J. Mod. Phys. D 6 563 (1997). 

[20] K. D. Patil and U. S. Thool, Int. J. Mod. Phys. D 14 873 (2005). 

[21] K. D. Patil and S. S. Zadc, Int. J. Mod. Phys. D 15 1359 (2006). 

[22] M. C. Bcnto, O. Bertolami and A. A. Sen, Phys. Lett. B 575 172 (2003); Phys. Rev. D 67 063003 (2003); 
Gen. Rev. Grav. 35 2063 (2003); D. Carturan and F. Finelh, Phys. Rev. D 68 103501 (2003); L. Amendola, 
F. Finclli, C. Burigana and D. Carturan, JCAP 07 005 (2003). 

[23] S. Pcrlmutter et al, Astrophys. J. 517 565 (1998); A. G. Riess et al, Astrophys. J. 116 109 (1998); P. de 
Bcrnardis et al. Nature 404 995 (2000); S. Hanany et al, Astrophys. J. 545 L5 (2000). 
[24] T. Padmanabhan, Phys. Reports 380 235 (2003). 



12 



[25] P. J. E. Peebles and B. Ratra, Rev. Mod. Phys. 75 559 (2003). 

[26] S. Chaplygin, Sci. Mem. Moscow Univ. Math, bf 21 1 (1904). 

[27] A. Kamenshchik, U. Moschella and V. Pasquier, Phys. Lett. B 511 265 (2001); V. Gorini, A. Kamenshchik, 
U. Moschella and V. Pasquier, gr-qc/MQ?>QQ2. 

[28] H. B. Benaoum, /iep-i/i/0205140; U. Debnath, A. Banerjee and S. Chakraborty, Class. Quantum Grav. 21 
5609 (2004). 

[29] M. C. Bento, O. Bcrtolami and A. A. Sen, Phys. Dev. D 66 043507 (2002); N. Bilic, G. B. Tupper and R. 
Viollier, Phys. Lett. B 535 17 (2001). 

[30] P. T. Silva and O. Bertolami, Astron. Astrophys. 599 829 (2003); A. Dev, D. Jain and J. S. Alcaniz, 
Astron. Astrophys. 417 847 (2004). 

[31] O. Bertolami and P. T. Silva, Mon. Not. R. Astron, Soc. 365 1149 (2006). 
[32] F. C. Santos, M. L. Bedran and V. Scares, Phys. Lett. B 646 215 (2007). 

[33] H. B. Sandvik, M. Tcgmark, M. Zaldraiaga and I. Waga, Phys. Rev. D 69 123524 (2004). 
[34] M. Marker, S. Q. dc Olivcira and I. Waga, Phys. Lett. B 555 1 (2003). 

[35] Z. K. Guo and Y. Z. Zhang, Phys. Lett. B 645 326 (2007), astro-ph/Q^Q&mi] astro-ph/QmQ7QQ; G. Sethi, 
S. K. Singh and P. Kumar, Int. J. Mod. Phys. D 15 1089 (2006). 

[36] K. Anguige and K.P.Tod, Annals Phys. 276 257 (1999); U. Mukhopadhyay and S. Ray, asiro-p/i/0510550. 
[37] S. W. Hawking and G. F. R. Ellis, (1973) The Large Scale Structure of Space-Time, Cambridge University 
Press. 

[38] W. Israel, Nuovo Cimento 44B 1 (1966). 

[39] U. Debnath, S. Nath and S. Chakraborty, Gen. Rel. Grav. 37 215 (2005). 



